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Abstract 

In this paper we study in details a system of two weakly coupled harmonic oscillators. 
This system may be viewed as a simple model for the interaction between a photon and a 
photodetector. We obtain exact solutions for the general case. We then compute approximate 
solutions for the case of a single photon (where one oscillator is initially in its first excited state) 
reaching a photodetector in its ground state (the other oscillator). The approximate solutions 
represent the state of both the photon and the photodetector after the interaction, which is 
not an eigenstate of the individual hamiltonians for each particle, and therefore the energies 
for each particle do not exist in the Copenhagen interpretation of Quantum Mechanics. We 
use the approximate solutions that we obtained to compute bohmian trajectories and to study 
the energy transfer between the two particles. We conclude that even using the bohmian view 
the energy of each individual particle is not well defined, as the nonlocal quantum potential 
is not negligible even after the coupling is turned off. 

1 Introduction 

The discussions about the incompleteness of the wavefunction to describe physical processes dates 
back to the beginning of quantum mechanics itself. This discussion is closely related to the 
possibility of describing quantum mechanical systems from an underlying realistic model. In 1952, 
David Bohm showed that such a realistic model was possible. However, Bohm's theory had the 
problem of being nonlocal QlEj. In 1963 John Bell showed that in order to obtain the same results 
predicted by quantum mechanics, any realistic theory would have to be nonlocal Bell's result 
and the failure of using the Copenhagen interpretation of quantum mechanics to some particular 
situations, as in for example Quantum Cosmology, lead to a raised interest in Bohm's interpretation 
and in nonlocal realistic theories (3. 

The subject of reallity and nonlocality has been an interest of Patrick Suppes for quite a while 
[T7| . in particular for the photon. In fact, one of the authors of this paper co-published with 
him a series of papers that layed down the foundational analysis of realistic and local model of 
photons that could explain the double slit experiment, the EPR experiment and other phenomena 
|18[ I19L I2H . The problem with the Suppes and de Barros model was that it did not have 
a consistent theory of photon-counting for single photons, and therefore could not explain the 
non-locality of single photons and the GHZ experiment, for example. 



In this paper we try to respond, within Bohm's model, the question: what is a photon? We 
do not follow the standard Bohmian interpretation for bosonic fields (as can be found in [TT]1. 
Instead, we use the simple interpretation that "a photon is what a photodetector detects". One may 
think of a photodetection as a transfer of energy from a quantized mode of the electromagnetic 
field (the photon) to an atom in its ground state (the quantum photodetector). Therefore, to 
study this photodetection we will focus on the process of transfer of energy from the photon to 
the photodetector. 

To study the exchange of energy in details, we have to choose between two different and 
simple models of a photo-detector: a photo-detector with discrete or continuous band p]. For 
the purpose of simplicity, we will choose the former. However, since we are only interested in the 
aspects of energy transfer between the two systems, we will make an even further simplification and 
consider that the photon and the detector will both be described by a single harmonic oscillator. 
Furthermore, during some time ATi„t, we will assume that a linear interaction exists between the 
two oscillators, and that this interaction is weak. This detection model is known as an indirect 
measurement and has been the subject of intense research lately as it is directly connected 
to quantum nondemolition experiments. As we will see, this "toy model" will allow us to capture 
some important features of the entanglement between the two systems. 

This paper is organized in the following way. In Section[2lwe will quickly review the interaction 
between two harmonic oscillators for the classical case. This will allow us to understand how the 
transfer of energy happens in such case. We then compute the exact solutions for the quantum 
mechanical system with interaction (Section I^J. In Section ^ we apply the results of Section to 
a specific case of exchange of a single quantum of energy and analyze its outcomes. In Sectional 
we use Bohm's theory to interpret the results obtained. The conclusions are in Sectional 



2 The Classical Case 

Before we go into the details of the quantum mechanical examples, let us begin by analyzing the 
classical system of two one-dimensional coupled harmonic oscillators with the same mass m, elastic 
constant k, and coupHng constant A, as shown in Figure Q The Hamiltonian for this system is 
given by 

^ = i (^i' + P2) + \k ((^1 + df + (^2 - df) + (Xi - X2 + 2df . (1) 

To simplify the equations of motion and eliminate the undesirable constant d we can make the 
canonical transformation 

x\ = A^i + d, 
X2 ~ X2 ~ d, 
PI = Pi, 

P2 = P2- 

With the new variables equation Q rewrites to 
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Figure 1: Identical harmonic oscillators coupled by a spring of constant A. 

H=^{pI+ pD + \k {^l + + {x, - X2f . (2) 
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The Hamiltonian equations of motion are 



Pi = -t; — = -kxi - \[xi - X2) , 
OXi 



Xi 



dH _ pi 
dpi m ' 



dH ^ 

P2 = -T; = -kX2 + X{Xi- X2) : 

0X2 



X2 



dH _ p2 
dp2 m ' 



yielding 



and 



m{xi+X2) = -A: (xi + X2) , (3) 



m{xi-X2) = - [k + 2X) {xi - X2) ■ (4) 
The general solutions to I^J and Q are 

\/2^+ =xi+X2= A cos ( y^t+^ ) > 



\/2C_ = xi - X2 = A' cos "4^* ^ ^' j ' 
and ^_ are the normal coordinates of the coupled harmonic oscillators) or, equivalently, 



xi = ^ cos I y + 6> j + Y cos \ \ r ' "" ^ + 0' 



Ik 


+ 2A 


m 


>k 


+ 2A 



We will assume that the two oscillators are initially at rest the first one in its equilibrium position 
(null initial energy, E\ =0), while the second one is dislocated from its equilibrium position by a 
distance D (initial energy given by E\ = {1/2) kD"^): 

ii(0) = i2(0) = 0, 

a:i(0) = 0, 

X2(0) = D. 

The integration constants then read 



-6' = 0, 
A = D, 
A' = -D. 



yielding 
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Xlit) 



D 



D 



COs(tjt) — COS {iij't)] 



X2{t) = — [cos (ujt) + COS (l^'t)] . 



(5) 
(6) 



where we defined iv = y/k/m and uj' = lu^/I + e, with e — 2X/k. Equations 10) and © can be 
written in the following suggestive way. 



xi{t) = — Z^sin 



X2{t) ~ D cos 



{u - Lu')t 



{uj - uj')t 



cos 



(w + iu')t 



{uj+uj')t 



(7) 



(8) 



We will now assume that the interaction constant A is weak when compared to the elastic 
constant fc, £ ^ 1. Then, we can expand uj' around e = 0, yielding 



k + 2X 



with 



Defining 



Suj — 



u}^l + e = ijj + 2Suj 



X 



2 2Vkrn 
uj' + u; 

— - — = UJ + OUJ, 



(9) 

(10) 

(11) 



the solutions can now be written as 



xi{t) = D sui{5uj t) sm[uj t], 
X2{t) = D cos{Sujt) cos[ujt], 



(12) 
(13) 



where the dependence on A of Eqs. lfT2)l and ifT^jl are present in Suj and uj through ((T(1|l and ((TT|l . 

The movement of both particles is periodic, with two characteristic frequencies Suj and uj. The 
frequencies Suj and u) are known as the normal modes of vibration, with uj being called the higher 
normal mode and Suj the lower normal mode. Both movements have period r = 27r/(I' and are 
modulated by a variable amplitude with much greater period given by r = 2tt/5uj. They are 
■jt/2 out of phase. We can compute the energy of the two particles, Ei = p\/2m + kxl/2 and 
E2 = pl/2m + kxl/2. They are 



2 
kD^ 



sin^ {Suj t) 
cos^ [Suj t) 



1 + 4— cos^(<Zit) 

UJ 

1 + 4— sm^(<Iit) 

UJ 



(14) 
(15) 



Due to the coupling, the particles exchange energy between themselves periodically, with period 
T = 2it/Suj. Each of the oscillators achieve its minimum energy value when the other have its 
maximum value. The maximum value of the energy can be a little bit bigger then kD^ 12. This 
may seem odd, but we must remember that the extra energy is due to the interaction energy 
A(a;i — X2Y /2 = ke{xi — a;2)^/4. It is easy to check that if we add this interaction energy to the 
sum El + E2 we obtain the total energy of the system 



kD'^ 



1 



Suj 



OiSuj% 



(16) 
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a value that is constant for the whole movement, as we should expect. For more details, see 
Refs.|8l[22]. where this system and generaHzations of it are analyzed with detail. Of course, as the 
Hamiltonian is time independent, energy is always conserved. 

It is also interesting to note that the total energy of the system depends on the coupling 
constant, as shown by Ijlfill . A quick analysis of the origin of the "extra" energy shows us that this 
happens because of the initial conditions chosen. The initial conditions from which we obtained 
ifTfil have the particle represented by X2 off its equilibrium position, whereas the other particle is 
at its equilibrium position, with both particles having zero kinetic energy. This initial condition 
obviously imply that the coupling spring, with elastic coefficient A, is also stretched from its 
equilibrium position, and therefore has nonzero potential energy at i = 0. If we use other initial 
conditions, the "extra" energy due to coupling does not appear. For example, we can choose both 
particles at an initial position where all spring have no potential energy (in our case, xi = X2 = 0) 
and one of the particles has some kinetic energy while the other particle has zero kinetic energy. 
With this set of initial conditions, the energy transfer from one particle to the other is the same 
as before, but no coupling energy is present in the total energy. 

3 Quantum Evolution: Exact Solutions 

Now we want to study the quantized version of the resonant spinless one-dimensional coupled 
harmonic oscillator presented in the previous Section. First we note that the total Hilbert space 
H, = Til ® 7^2 is spanned by Tii and 7^2, the Hilbert spaces for particles 1 and 2, respectively. For 
example, the two canonical variables describing particle 1 are 

with 

and are therefore represented as 

where I2 S 7^2 is the identity operator. In this way, the Hamiltonian operator for particle 1, is 

written as 22 

Hi = ^ (Pi ®l] + ifc fXi (g) i + fii ® . 
2m \ / 2 V / 

For shortness of notation, we will drop the tensor product and keep in mind that operators 
regarding particle 1 act on Tii whereas operators regarding particle 2 act on H2 . 

With the simplified notation, the total quantum Hamiltonian operator for the two oscillators 
plus the interaction term is 

H = H1+H2+H1 

We can now make the following change of variables, similar to the classical case: 

xi = -^1 + d, 
X2 = X2 - d, 
Pi = A, 
P2 = Pi- 

This change of variables obviously keeps the commutation relations between momenta and posi- 
tions. Hence, in the coordinate representation we have the Hamiltonian operator 
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In analogy to the classical case, we work with the normal coordinates 



1 , 

T2 " • 



(19) 
(20) 



This change of variables has Jacobian one, and does not change the normalization of wavefunctions. 
With the normal coordinates, the Hamiltonian is 



H = - 



2m yd(l 



and is now separable, i.e. 
where 



H = H+ 
2m d^l 



(21) 
(22) 

(23) 
(24) 



Equations l(23|l and ll23|) are the well known Hamiltonians for one-dimensional uncoupled harmonic 
oscillator with frequencies ^Jk/m and \J{k + 2\) /m, respectively. 
The Schroedinger equation for the system is 



(25) 



To solve l|25ll we need to find its eigenfunctions and eigenvalues since they form a basis for the 
Hilbert space. The general solution can be written as a superposition of the eigenfunctions. Hence, 
we need to find the solutions to the time independent Schroedinger equation 



i?V'<''(e+,e-) = fiV'<"(e+,e-), 



(26) 



where I is an index (perhaps a collective index for both oscillators) for the energy to be determined. 
Since H is separable, we can write ll26|l as two independent eigenvalue equations 



i7+0^)(e+) = i?„0?^(e+) 



and 

where we define 
and 



V'(')(e+,e-) = 0V"^(e+)0L"'\e-), 



in'). 



(27) 

(28) 
(29) 



Clearly, I is an index that depends on both n and n', and for that reason we will write )(^_|_,^_) 
instead of ^''-'HC+j ?-)• The eigenfunctions of lf27jl and lf28l are well known to be 



1/4 



1/2 



7r?i22»(n!)2 



H„ 



exp 



2h 



(30) 



6L"')(e-) = 



{k + 2A) 



1/4 



7r?i22»'(n'!)2 



Hn' 



^m (k + 2A) 



1/2 






exp 



Vm(fc + 2A)C2' 



2h 



(31) 
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and its corresponding eigenvalues are 



En 



and 



E', 



k + 2X 



(32) 



(33) 



where Hn are the Hermite polynomials of order n P] . 

The solution to the time dependent Schroedinger equation H25II is obtained applying the time 
evolution operator 

U{t, to) = cxp (-iH{t - to)/h] 



on tjj{^+,^-,to). Since C-) = 0? ^(?+)0-"''(C-) fomi a basis for H, we have 

OO 



OO pOO 



^V"' (e+)</'^''^ (e- )V'(e+, , io) de+de- , 



— OO ^ — OO 



and we used the reality of \^-) in the expression for Cn,n' ■ Then, 

^{i+,i-,t) = ?7(t,to)^(^+,e-,0) 

OO 

n.n'— 

OO 



(34) 



= and we assumed, 



where exp = exp ^— ti/+</7i^ cxp ^— zif-i/?!^ since Hj^,H^ 

for simplicity, that to = 0. 

We can now finally go back to the original coordinate system xi and X2, and the explicit form 
for the general solution in this coordinate system is 



%j}{xi,X2,t) = 



1/4 



muj\ 
~2h) 



1/2 



22"(7i!)2; y22"'(n'!)2 

„ /\ 1/2 



1/4 



(xi + X2) 



2h 



(xi - X2) 



exp 
exp < —I 



r m r , ,2 1 , 1 

lihr +LU(Xl-X2)^J 



+ -]uj+ \n' + 



(35) 



where we defined, as before, uj ~ ^kjm and tj' = ^Jk + 2X/rn . The wavefunction l(35]l thus 
describe spinless one-dimensional coupled harmonic oscillators with no approximation. 
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4 A Simple Example 



We saw in the classical case that two coupled oscillators can transfer energy to each other. This 
was clear with the example where at t = one oscillator had zero mechanical energy while the 
other one had nonzero potential energy. As time passes, the mechanical energy of the former is 
transferred to the latter. It is interesting to study the quantum mechanical analogue to this case, 
i.e., when one quantum oscillator is in an excited state and the other is in its fundamental state. 
We will consider as the initial wavefunction the following 




X2 exp 



2h 



The wavefunction is an eigenstate of the Hamiltonian 



H = H, 



Ho 



(36) 



(37) 



without the interaction term Hj. Clearly, ^(2^1, 2:^2, 0) is separable, i.e., it is possible to write 
tlj{xi,X2,0) = (pi{xi,0)(p2{x2,0). Since Hi {H2) acts only in (pi{xi,0) {(p2{x2,0)), the state 
ifj{xi,X2,0) represents a system where the particle described by xi is in the ground state and 
the particle described by X2 is in the first excited state. So, we can think of our example as the 
following. We have initially a system of two harmonic oscillators, one in the ground state and 
the other in the first excited state. After t = we suddenly turn on a interaction between the 
two oscillators, and as a consequence we expect to have a "transfer of energy" from one oscillator 
to the other, as it happens in the classical case. We will now proceed to analyze in details this 
example. 



4.1 Approximate Solution 

To use equation (|35|l we need to find the coefficients C„„' . It is straightforward to compute the 
coefficients from l(3ijl by just using the orthogonal properties of the Hermite polynomials and by 
rewriting Il3(i|l in the normal coordinates, yielding 



Cn,n' — 



1/4 



22"(n!)2 

-rOn' ,2j(l\. . 



22"'(n'!)2 



1/4 



(2j 



D! 



(cj + t^') j! 



2w' 



7'5n',2j + l'5o,i 



(38) 



where i5y is Kroenecker's delta. 

It is interesting to note that there exists infinite terms of Cn,n' that are different from zero. 
Therefore, if we write down the expression for the time evolution of the wavefunction after the 
interaction we obtain an expression with an infinite number of terms. However, a close look at 
the Gn,n' coefficients may shed fight on how to deal with this problem. First we see from l(38|l 
that only the terms Co,,i' and Ci^„' are nonzero. If we compute the ratio between two consecutive 
nonzero terms, i.e, Co,Ti'+2/C'o,n' and Ci_„'+2/C'i,ti' we obtain 



Co n' 



+2 



Co,i 



Cln' 



+2 



C'l.n' 






^2) 


{n' - 


M)' 


\{n'- 


f 1) 



(n' + 2)' 



(39) 



(40) 



We note that both ratios l(39jl and ijifljl are proportional to 



Then, if the coupling constant 



A is small compared to k (weak coupling) we can make an expansion of l(39| and around A = 



8 



and obtain, up to first order in A, that 



A 



0(A2 



We conclude that if A is small compared to fc, as we increase the value of n', the coefficients Cn,n' 
become less important. Therefore, it is justifiable to keep only a finite amount of terms in the 
expression for ■ip{xi, X2,t) for small A. In our example, we will keep only terms up to first order in 
A. 

Since we will be working with A small, it is convenient now to introduce the following parameters 
already used in the classical case 



6lu 



OJ 



2 

Uj' + LU 



Then, if A is small. 



and 



Suj 



< 1. 



Keeping only terms up to first order in we have 



7i,n'—0 



(41) 



where 



Cio 
Coi 

Cl2 
Co3 



V2 
2 ' 

_V2 

2 ' 
1 Slu 

a/3 6lij 



(42) 

(43) 
(44) 

(45) 



We are finally in a position to write, up to first order, the time dependent wavefunction for the 
coupled harmonic oscillators. From lUTjl and (Il2)l - (|l5ll it is straightforward to obtain 



, , I 1 muj r muj r , 21! r r,- i 

tp{xi,X2,t) = d——exp^-—[x^+X2\jexp{~i2ujt} 



2i { xi + X2 



1 muj 2 

2 """"1 



^ I sin (Sujt) + 2 [ X2 + xi 



1 mu) 2 

2 """"^ 



Suj 



cos {5u)t) 



1 5uj , . \mUJ \2 , r m- ^ \ ^ 

+ -— (a:i+X2) -^{Xi-X2) -1 exp{-i(2a; + dw)i| 

2 UJ in 



1 Suj r ttluj 

---^(xi-xa) ^{xi~X2f cyi^{-i{2u) + i5uj)t]\ +0{5ijj'^). (46) 

2 UJ \- Ti J 



The wavefunction (I46|l determines the evolution of the system. We will now proceed to analyze 
the system using 
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4.2 Marginal Probabilities 

From H4(ill we compute the joint probability density for xi and X2 as a function of t. The joint 
density is simply 

Pixi,X2,t) = \^{xi,X2,t)\'^, 

and keeping terms up to first order in Suj we have 

exp|-^ [2:? +2:2]} X 

1 mu) 

2 ~ T 
1 miu 



P{xi,X2,t) 



1 / muj \ ^ 



2tt \ h / 



4 yx2 + 2xiX2 

+4:(xl + 2xiX2 



-X2 



-x\ 



^ ) cos^ ((5a;t) 
— ) sin^ {6u}t) 



5uj , 
+2x2 — (2:1 

OJ 


+ X2) 


'muj 
-~ 


{xi 


- X2) 


-2X2— \X\ 

Ijj 


- X2) 


'muj 
-~ 


{xi 


- X2) 


—2x\ — [x\ 

UJ 


+ X2) 


'muj 
-~ 


[xi 


- X2) 


Sll> 

+2xi — [xi 

UJ 


" X2) 


'muj 


{xi 


- X2) 



It is interesting to see how the marginal probability distributions for xi and X2 behave. Let us 
recall that the marginals are defined as 



cos {6ujt) cos {{2uj + iSuj) t] 
sin {5ujt) sin {{2uj + 5uj) t} 
sin {5ujt) sin { {2uj + 2,6uj) t] \ (47) 



P{xut) 



P{xi,X2, t) dx2, 



and 



/oo 
P{xi,X2,t) dxi. 
'OO 



(48) 
(49) 



Therefore, P{xi,t) dxi represents the probability of measuring the position of particle 1 in the 
interval {xi^xi + dxi) independently of particle 2. The interpretation for P{x2,t) is similar. 

From Ij47ll . Ij48|l . and II49|I it is tedious but straightforward to compute (once again up to first 
order in Suj) such quantities, which read 



P(xi,t) 



muj 



( muj 2I 
exp|-_xi| 



2 , - ,, 2muj 2 . 2 / r ,\ 
cos (Sujt) H T—xf sm [Sujt) 



Suj 

UJ 

muj 



1 mui 2 



muj 



(3 cos ((2w + 3Suj) t) — cos {{2uj + Suj) t)) cos (Sujt) 
(sin {{2uj + 3Suj) t) - sin {{2uj + Suj) t)) sin {Sujt) 



(50) 



and 



P{x2,t) 



( muj 2} 
exp|-_x2| 



ttTl 



sm (Sujt) H — X2 cos [oujt) 



Suj 

UJ 

muj 



1 muj 2 
4 ~ ^""2 

3 muj 



2 h 



(3 sin ((2w + 'SSuj) t) + sin {{2uj + Suj) t)) sin {Sujt) 
(cos {{2u) + 3(5tj) t) + cos {{2uj + 5w) t)) cos (Jcjt) 



,(51) 
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We can compute the values of the marginals and (|5Tll at i = and find that, after making 
sure that we use uj as the frequency instead of w, and keeping only terms up to first order in Jw/w, 
such marginals indeed represent the ones for the ground state HO and the first excited state HO, 
as one should expect. 

To better grasp the behavior of 15011 and Ij51|l . let us plot them as a function of time. Before 
plotting, we need to choose the appropriate values for the constants in the equations. If our 
system is in atomic scale, it is not reasonable, from a computational point of view, to use the MKS 
system. So, we will measure time in femtoseconds (Ifs = 10~^^ s) and distance in Angstroms (1 
A= 10"^" m). If we say that the particles in the oscillators are electrons, then m = Ime, where 
me is the mass of the electron, then we have 

n = 10me-A^-fs"\ 

and 

fc = 1 me • fs~^, 

and, for the harmonic oscillator. 




Figure 2: Graphs for the marginal probabilities of xi and X2 as a function of time. In these graphs 
we used m = 1 me, u) ~ 1 fs~^ and 5uj/u} = 1/10. The scale for time is fs and the scale for distance 
is A. 

((Ax)2) = 

The behavior of the probability density for particles 1 and 2 are found in Figure[2| The time interval 
chosen for the time axis in the graphs was A< = tt/cJoj as this is the value where cos {Sujt) = — 1, 
which is an extreme in the behavior of the marginal densities. Looking at the graphs we see that 
particle 1 starts with a marginal density that is mainly a Gaussian function, whereas particle 2 
starts from the product of a;| times a Gaussian. This is because particle 1 is at the ground state 
and particle 2 is at the first excited state aX t = 0. However, as time passes there is a swap in the 
roles of particle 1 and 2, in the sense that at t = tt/Suj the marginal density for particle 1 resembles 
that of particle 2 for i = and vice versa. This is of course due to the interaction between the 
two particles. We may think of those densities as showing that, at t = tt/Suj (more generally 
when t = {2n + 1) tt/Suj) particle 1 is no longer in the ground state, but in the first excited state, 
whereas particle 2 is in the ground state. 

4.3 Energy Expectations 

The densities above suggest that there is an energy transfer from one particle to the other. To 
see that this is the case, let us compute the energy values for each particle. First we should note 
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that the system is not in an eigenstate of the Hamiltonian, as we started from a superposition of 
different energy states. We define the energy or particle 1 as 



El = (Fi), 

the energy of particle 2 as 

E2 = {H2), 

and the total energy as the sum of the two energies plus the interaction energy 

Et^Ei+E2 + {Hi). 
In coordinate representation we have that 

/>oo />oo 

El = / dxidx2il^{xi,X2,t)*Hiipixi,X2,t) 



— 00 — 00 



00 />oo 



dxidX2 "0(2^1: ^2, t)* 



00 J —00 



h^d^ 1 2 

'2mdxl ^ 2 \ 



i>{xi,X2,t), 



and computing this term we obtain, up to second order in Suj/lo, 

El = hu! ( ^ + sin^ {Sujt)] {1- Suj/lu) 



hLu{^ + sin^ {Sujt) ] . (52) 



Similarly, for E2 we have 



E2 ^ hu! + cos'^ (Suit) ] {l-6uj/oj). 



= ^it^ Q + cos^ (SLut)^ . (53) 

If we compare the quantum energies (|52ll and (l53)l to the classical expressions ifTHl and (tTfijl the 
resemblance is striking. They are practically the same for Suj/Oj <^ 1, except for a zero energy 
factor of ^hOj present in the quantum mechanical case. In fact, the same conclusions can now be 
drawn from II52|I and (l5;-i|l , i.e., that due to the coupling, the particles exchange energy between 
themselves periodically, with period r = 2Tr/Suj. Each of the oscillators achieve its minimum 
energy value when the other have its maximum value. For the interaction energy we compute 

(Hi) = 2h5uj. (54) 

Then, it is easy to compute the total mean energy 

Et = El -\- E2 + {Hi) 

= 2hLU 

= 2huj + 2n5uj. 

This is once again in agreement with the classical case seen above, in the sense that the total 
energy is the sum of the energy of each oscillator (keeping into account the nonclassical zero point 
energy) without the interaction term plus an interaction term 2fi5uj. 

We just saw that the state we used had a term in the total energy 2Ti5uj that was due to the 
coupling between the two oscillators. However, if we remember the classical case of Section 2, with 
different initial conditions — e.g. xi = 0, X2 = 0, ii = w, ±2 = 0, at t = — no interaction term is 
present in the total energy. What about the quantum case? Do we always have an interaction term 
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present, as in (Ifiljl ? A short computation shows that for any initial state that is a combination of 
Fock states for the two HO of the form 



IV') = \ni) ® |n2), 

where |ni) and \n2) are eigenstates of two uncoupled HO, the value of (Hj)^ (the interaction 
term) is different from zero. 

The question remains as to whether it is possible to find an initial state that has an interaction 
term that is zero. A good guess would be to take both HO in a coherent state at t = 0, since it is 
a state that has many of the characteristics of a classical system It is easy to show that it is 
indeed true that for the state 

where 

\a) = e 



n— 



and similar for the expected value of the interaction energy at i = is zero if a and /? have 
an appropriate phase relation. It is left up to the reader to find out this phase relation and a 
set of initial conditions for a classical system which reproduces the expectations in the quantum 
mechanical case. 



5 The Bohmian Interpretation 

Before we analyze the transfer of energy from a Bohmian point of view, let us quickly review 
Bohm's interpretation of quantum mechanics. Let us begin with the causal interpretation for the 
case of the Schrodinger equation describing a single particle. In the coordinate representation, for 
a non-relativistic particle with Hamiltonian H = fp' /2m + V(x), the Schrodinger equation is 



*(a;,t). (55) 



We can transform this differential equation over a complex field into a pair of coupled differential 
equations over real fields. We do that by writing = Rexp{iS/h), where R and S are real 
functions, and substituting it into l|55ll . We obtain the following equations. 

f +<Z^)! + V-f^ = 0, (56) 

ot 2m 2m R 

^+y-{R' — )^0. (57) 
ot m 

The usual probabilistic interpretation, i.e. the Copenhagen interpretation, understands equation 
ll57jl as a continuity equation for the probability density R^ for finding the particle at position x 
and time t. All physical information about the system is contained in i?^, and the total phase S 
of the wave function is completely irrelevant. In this interpretation, nothing is said about S and 
its evolution equation . 

However, examining equation (j57ll . we can see that VS/m may be interpreted as a velocity field, 
suggesting the identification p ~ \/S. Hence, we can look to equation lf5fijl as a Hamilton-Jacobi 
equation for the particle with the extra potential term 

2m R ' 

where Q is the so called quantum potential. Thus, since Bohm's interpretation identifies p with 
VS*, from the differential equation p = mi = V5 we may compute its solutions and obtain the 
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trajectory of the quantum particle. Therefore, in Bohm's interpretation both momentum and 
position are quantities that are ontologically well defined. 

For our case of two coupled-HO, the configuration space has two variables, xi and X2, repre- 
senting the positions of particles 1 and 2, respectively. For two particles, the nonlocality of Bohm's 
interpretation becomes evident as the Schrodinger equation becomes 



di 



2m 



2to2 



V{xi,X2) 



'^{xi,X2,t), 



(58) 



where V| is the laplacian operator with respect to the coordinates of particle i. If we follow the 
same transformation as before, we can obtain the following equations. 



95 
dt 

die 
dt 



2m2 



2mi 
-Vi • 



V - 



2m R 













V "12 / 



= 0, 



0. 



(59) 
(60) 



The nonlocality comes from the fact that, even if the potential V{xi,X2) is local, it is possible 
that the quantum potential given by 



Q = - 



vlR 

2mi R 2m2 R 



VlR 



is nonlocal, depending on the form of R. This characteristic is necessary, as proved by Bell, if 
Bohm's theory is to recover all quantum mechanical predictions. 

Using (|46ll it is straightforward to compute the phase S{xi,X2.,t) from the expression 

^{xi,X2,t) ~ ^{xi,X2,t)* 



S{xi, X2,t) = — ?i arctan 
After some long and tedious algebra we obtain 

S{xi,X2Tt) = —h arctan 



^{xi,X2,t) + '${xi,X2,t)* 

SA{xi,X2t) ^ 
SBixi,X2t) 



where 



SAixi,X2,t) = A COS {2Lut) ^(^xl sin {6ujt)^ — X2 COS (SLot)"^^ sin {2bjt) 
+X1X2 sin {Sujt) cos (Sujt)} , 



and 



SB{xi,X2,t) 



{x2 cos{2ujt) cos{Sujt) + xi sin(2a;t) sin{dujt)) , 



where we keept all terms in {Slu/uj) t but we neglected terms in Suj/Cj. 

From S(xi,X2t) we obtain the differential equation that describes the trajectories of particles 
xi and X2 as 

Axi 1 dS{xi,X2t) h X 2 cos (Sujt) sin {6 ujt) 

(61) 



and 



dt 

dX2 

"dF 



71 dxi m x\ sin? {Scot) + x\ cos^ {6ujt) 

1 dS{x\^X2() h X I cos (dujt) sin {S Lot) 



dx2 



(62) 



m x\ sin^ (Stot) + x\ cos^ [bu^t] ' 

We can see that the trajectories follow a set of differential equations that are coupled and nonhnear. 
It is interesting to notice that if iJo; = we recover the standard Bohmian result that in the case 
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of no interaction each HO is in an eigenstate and therefore both particles are at rest. However, if 
5u ^ 0, we obtain at once that, after the change of variables 




t' = — t (63) 

(64) 
(65) 

the diflerential equations (I61|l and Ij62|l are form invariant with respect to a change in the coupling 
constant from 5uj to Soj' . This invariance is illustrated in Figures El and ^ where typical Bohmian 
trajectories were computed for both particles. The solutions shown in Figures El and 21 were 
obtained numerically using a 7th-8th-order continuous Runge-Kuta method. 
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Figure 3: Bohmian trajectories for two CHO. The trajectories correspond to w = Ifs"'^, Slo/uj = 
0.01, a;i(0) = 0, and X2(0) = —1. The solid line represents the trajectory of xi{t) whereas the 
dashed line represents that of X2{t)- The scale for the ordinates is in A and the time scale is in fs. 



It is important to compute, in Bohmian theory, the quantum potential Q defined as 



where 

^ 1 aVP(.Tl,X2,f) 

2to ^P(a;i,X2,t) dx\ 
and 



2to y^P{xi,X2,t) dxl 

It is straightforward to compute 

Qi = fiu) — —milp'x\ 

^ 1 ^_ sin^ {5ujt) — X2 cos^ {Sujt) 
2 x\ SVC? {Scot) + x\ cos^ {6ujt) 
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Figure 4: Bohmian trajectories for two CHO. The trajectories correspond to w = lfs~^, (5w/w = 
0.005, a;i(0) = 0, and a;2(0) = — \/2. The soHd Hne represents the trajectory of xi{t) whereas the 
dashed Hne represents that of X2{t). The scale for the ordinates is in A and the time scale is in 
fs. We can observe that the trajectories are identical to the ones shown in the previous Figure, 
except for the coordinate scales, a result consistent with equations llfi3)) - (|65ll . 



1 fi^ x\ cos^ {Suot) sin'^ {5ujt) 



2 m 



'x\ siv? {8ujt) + x^ cos^ {5ujt)y 



and 



^ -2_2 



Q2 = hu) — -mu! X 



1 ^_X2 COS^ (Scdt) — x\ SVC? (Sbjt) 

2 X2 cos^ (Sut) + x\ SIT? (Sujt) 
1 ?? x\ SVC? (Sojt) cos^ (Sujt) 



which yields 



(66) 



(67) 



Q{xi,X2,t) = 2huj ~ -mu!^ (^x\ + X2) 

1 ?? {xl + x'^) sir? (Sujt) cos^ (Suit) 

We are now in a position to compute the total bohmian energy for each one of the particles, 



El = Ki 



Vi 



1 ^ sin^ (6ujt) — x'i cof? (6ujt) 

= hQ + -hu)-^ — 2^ — ^ i — 

2 xf sin {Scut) + X2 cos^ {Scot) 



E2 = K2 + V2 + Q2 
1 ..2 

= huj — huj 
2 



^^- xl sir? (Sut) — X2 cos^ (Sujt) 
2 xl sir? (Sujt) + x\ cos^ {6ujt) ' 
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where Ki = \m is the kinetic energy of particle i (obtained from the guidance equations 

iffiTll and ijESl) and Vi is the potential for particle i (neglecting terms in 5uj/uj). 

The total energy for the system is just the sum of the individiual energies, yielding 

Et = Ei+E2^ 2nu!, 

the same value as the expected energy of the system. 

6 Conclusions and Final Remarks 

We see that the expressions obtained for Ei and E2 involve an interaction term that makes it 
impossible to distinguish what part of the energy belongs to the particle xi and what part belongs 
to the particle X2, except for some particular values of t. In the Copenhagen interpretation of QM 
it does not make any sense to talk about the energy of each oscillator for all t, as the oscillators 
are in a quantum superposition and are not in an eigenstate of its hamiltonian operator. In Bohm, 
it will not make any sense to talk about the energy of each oscillator for all t, since the quantum 
potential creates an interaction between the two oscillators that is of the same order of the other 
terms in the hamiltonian. Therefore it does not make any sense in the bohmian theory to say that 
the energy of the photon was transfered to the photodetector (except for very special values of t) . 

However, the bohmian interpretation gives an onthological explanation for the indefiniteness of 
the energy of each particle. Even with the interaction turned off, there is still a quantum nonlocal 
interaction between the oscillators given by the quantum potential and, in fact, one oscillator 
is not isolated from the other. This indicates that a real measurement has not yet ocurred. It 
seems to us that in order for a measurement to take place, a more elaborated description of the 
photodetection process involving a thermal bath or a macroscopic description must be used. In 
such case, we expect that the quantum potential will vanish and no further nonlocal interaction 
will be present after the measurement. 
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